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The baryon-number-two {B — 2) solution based on the SU(2) chiral quark soliton model (xQSM) 
is solved numerically, including fully the sea quark degrees of freedom. We confirm that the axially 
symmetric meson configurations yield the energy minimum for the B = 2 state in the xQSM 
when taking into account quark dynamics. Due to the axially symmetric meson fields, six valence 
quarks occupy the lowest energy level, consistent with the Pauli exclusion principle. The minimal- 
energy of the entire system is obtained within the framework of a symmetric ansatz. The fermion 
determinant with axially symmetric meson fields is calculated by diagonalizing the corresponding 
Dirac hamiltonian in a non-perterbative way, using a cylindrical Dirac basis. The baryon number 
density, calculated with quark fields corresponding to a soliton, is toroidal in shape. We also 
calculate the mean radius of the toroid from the quark fields. These results are closely related to 
Skyrme model calculations based upon pion degrees of freedom. Our model calculations clarify the 
underlying dynamical structure of the baryons at the quark level. 

PACS numbers: 13.75.Cs, 12.39.Fe, 12.39.Ki 


Although QCD is generally accepted as the underlying 
theory of the strong interaction, most low- and medium- 
energy nuclear phenomenology may be successfully de¬ 
scribed in terms of the hadronic degrees of freedom. In 
the case of the deuteron, the simplest nuclear system, 
similar approaches exist Recent investigations for 

deuteron photodisintegration and deep inelastic scatter¬ 
ing of leptons by nuclei suggest the necessity of including 
quark degrees of freedom^, ^ For this reason 

it was suggested that nuclear theory should be reformu¬ 
lated, taking into account the underlying theory. How¬ 
ever, QCD has difficulty in describing low- and medium- 
energy nuclear phenomenology because the coupling con¬ 
stants become extremely large at these energy scales, and 
it is desirable to formulate an effective, tractable theory 
for the strong interaction. 

The most important feature of QCD in the low energy 
region is chiral symmetry and its spontaneous break¬ 
down, or the appearance of Goldstone bosons. A sim¬ 
ple quark model that incorporates the above features is 
the chiral quark soliton model(xQSM)j|, pi] p^, 
which is charactrized by the following formulae: 


Z = 


'Dn'Dif'Dtp^ exp 




( 1 ) 


with , 

For the case of baryon number one {B = 1), the model 
was solved numerically in the Hartree approximation 
with a hedgehog ansatz: tv{x) = rF(r )&00E10- 
The profile function F{r) varies between the topologi¬ 
cal boundary conditions F{Q) = —tt and F{oo) = 0. A 
soliton solution with three valence quarks was obtained. 
The solution was identified as the nucleon (and A) after 
projecting onto good spin-isospin states. 

In the same way, the hedgehog ansatz in which the pro¬ 
file function varies from F{0) = —27r to F{oo) = 0 gives 


the baryon number two{B = 2). But the minimal en¬ 
ergy of this configuration was as large as three times the 
B = 1 mass[^. Therefore, the ansatz yields no “bound 
state” of the B = 2 system. 

On the other hand, in the Skyrme model the minimal 
energy configuration for B = 2 has axial symmetry |l4| 
H, 0 0 |i§. If we choose the z axis as the sym¬ 

metry axis, meson fields U with a winding number m are 
given in Manton’s conjecture|Q, 

U = cos F{p, z)ir ■ fismF{p, z), (2) 


where 

n = (sin 0(p, z) cosrrvp, sin 0(p, z) sinrnip, cos 0(p, z)). 

(3) 

The functions F{p,z), 0(p, z), the profile functions, 
are determined by minimizing the classical mass derived 
from the Skyrme lagrangian. The calculation was done 
by Braaten and Carson |19| . They obtained toroidal con¬ 
figurations which were classically stable. The ground 
state of this solution had the quantum numbers of the 
deuteron; however, the calculated static properties were 
not always consistent with the deuteron |p^. 

In this paper, we investigate the B = 2 soliton with 
axially symmetric meson fields in the xQSM. 

First, we introduce the one particle Dirac hamiltonian 
F[{IF'^) and a complete set of single-quark states as the 
eigenstates of this hamiltonian, given by 

= f3{idt - H{U ^^)), 


H{U^^) = -ia ■ V -k , 

(4) 


(5) 
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Using these eigenvalues, the total energy of the system 
can be written as 

^static [U] = NcEi^^[U]+N,Ei^^[U] 

+Efield[U] — Efield[U = 1] , ( 6 ) 


where 




0 -^0 ’ 


Ef^eld = ^A/;I + exp - (^) ^ 


with 


Nu = - 


•\/ 4 ^ \2 


(f) 


(7) 

C8) 

(9) 


The and Efield stand for the valence quark contri¬ 
bution to the energy for ith baryon and the sea quark 
contribution to the total energy, respectively. Here, no is 
the occupation number of the valence quark; that is, no 
is 0 or 1 . Efield is evaluated by the familiar proper-time 
reguralization scheme[^. A is the cutoff parameter. 

The Dirac hamiltonian E[ with axially symmetric me¬ 
son fields U commutes with the operator = L'i + 
For m = 2, H also commutes with the 
T 3 . Due to the symmetries of the me¬ 
son field configuration, above operator /3 • T 3 works as 
the parity operator. (For m = 1, the parity operator 
is given by a conventional form V = f3. ) L 3 , as and 
T 3 are the third-component of the orbital angular mo¬ 
mentum, the spin angular momentum, and the isospin 
operator of the quark, respectively. K 3 is often called 
the third-component of the grand spin operator. Conse¬ 
quently, the eigenstates of H are specified by the mag¬ 
nitude of K 3 and the parity tt = ±. As L 3 is integer 
and as and T 3 are ±1, so the possible values of K 3 are 
±|, ±|, ±1 • • • for m = 2. H also commutes with the 
“time-reversal” operator T = * 7173 - 1 x 1 T 3 C, where C is the 
charge conjugation operator. By virtue of this invariant, 
we see that the states of +K 3 and —K 3 are degenerate in 
energy!^, According to the Kahana and Ripka[p^, 


50-3 + ^rriTs 
operator V = (3 


we begin with investigating the spectrum of quark orbits 
as a function of the “soliton size” X (see Fig. 1). We find 
that only the Ks = states dive into the negative- 
energy region as X increases. Therefore one conclude 
that the lowest-lying axially symmetric B — 2 configura¬ 
tion is obtained by putting three valence quarks each in 
the first two positive energy states K 3 = . In that 

case, one immediately find that six valence quarks are all 
degenerate in energy. As a result, in our B — 2 system 
each baryon has equal classical mass. The degeneracy of 
the baryon in our system is a distinct feature of choosing 
axial symmetry as the symmetry of the meson fields. 

On the other hand, if we adopt the hedgehog ansatz 
for the B = 2 system, one find that the resultant hamil¬ 
tonian commutes with the grand spin operator K and 


the parity operator V. The magnitude of the grand spin 
operator K has the values of 0, 1, 2, 3- • •, then if we 
first put three valence quarks on the state AT = 0+, 
the next three quarks must be placed in higher energy 
states. If the second strong level AT = 0“ is occupied by 
next three valence quarks, the total energy is about 4260 
MeVQ. If we do not restrict the problem to “Skyrme- 
like configuration” , one can place the quarks into 
AT = 1+ or AT = 1“. Unfortunately, there are no thor¬ 
ough analysis which adopt these configurations within 
the framework of the yQSM. In the u-model calculation, 
their total energy is about 2620 MeV for the K = 0+, 1+ 
configuration ||24[ . In any case, the resulting total energy 
is much larger than two times of the mass with isolated 
baryon. As a result, we confirm that the energy of the 
axially symmetric K 3 = ± 5 configuration is lower than 
that of all possible B = 2 hedgehog configuration. Fi¬ 
nally we conclude that the lowest-lying B = 2 state has 
axial symmetry, and is obtained by putting six valence 
quarks in the degenerated states Ks = 

The meson field configuration that minimizes the to¬ 
tal energy Estatic\U] is determined by the following ex¬ 
tremum conditions: 


5E{p, z) 


Estatic\U] — 0 , 


5Q{p, z) 


Estatic\U\ — 0. (10) 


By using the explicit form of Estatic\U], this yields the 
following equations of motion for the profile functions 
A(p,z) and 0 (p,z), 


Ai2(p,z)cos0(p,z) = i?3(p, ^)sm0(p,z), (11) 

S{p,z) cos F{p,z) = P{p, z) sin E{p,z), (12) 


and 


P{p, z) = i?i 2 (p, z) sin 0(p, z) + Rs{p, z) cos 0(p, z ), 

(13) 

where 


Ri 2 {p,z) = 2Ri2v{p,z) + Ri 2 s[p,z), (14) 

Rs{p,z) = 2Rsy{p,z) + Rso{p,z), (15) 

P{p^z) = 2P^{p,z) + Pq{p^z) , (16) 


where subscripts v and 0 denote the contributions from 
the valence quarks and the sea quarks, respectively. The 
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FIG. 1: Spectrum of the quark orbits are illustrated as a function of the “soliton size” X. Profile functions are given by 
F{p, z) = —TT + 7r(p^ + /X, and 0(p, z) = tan“^(p/z). The orbits are labeled by . 


explicit forms of R 12 , R 3 and S are 
Ri2vip,z) = no 

X / dip$o{p,P,z)i^ 5 {TiCOS(fi + T 2 Smip)(j)o{p,<fi,z) , 

(17) 

Ri2o{P, z) = E Mvsign{E,j) 

V 

X / (i<p(^,y(p,(p, z)z75(ri costp + T2 sin(p)(/)y(p, (p, z), 

( 18 ) 

R3v{p,z) = no j dif^o{p,P,z)n5T34‘o{p,P,z), ( 19 ) 
R3o{p,z) = y^JV„sign{E^) 

V 

X J dip^^{p,(p,z)i'j 5 T 3 (j)^{p,(f,z) , (20) 

Sy{p,z)=no j dip^o{p,P,z)(l}o{p,p,z), ( 21 ) 

Soip,z)= y^Af^signjE^) 

U 

X j dip^u{p,ip,z)(j)^{p,ip,z), ( 22 ) 

In order to evaluate Eqs. ( 11 ) and ( 12 ) numerically, 
the following procedures were employed. Firstly, we 


start from the initial functions Fo{p, z) and 0o(p, z) that 
satisfy the boundary conditions given by Braaten and 
Carson [p^: 

F{p,z)^0 as p^ + ^ 00 , (23) 

F(0,0) = -7r, 0(O,z) = |“’ (24) 

[_7r, z < U 

We solve the one particle Dirac equation using the above 
functions Fo(p, z) , 0o(p, z). Secondly, i?i2(p, z), i?3(p, z) 
and S'(p, z) are calculated from the resultant eigenvalues 
and eigenfunctions; Q{p,z) is given by Eq. (11). Thirdly, 
the function F{p, z) is obtained on the basis of Eqs. (12) 
and (13). Then, new iterates of F{p,z) and Q{p,z) are 
obtained by resolving the Dirac equation. This procedure 
is continued until self-consitency is attained. 

Before reporting our results, we provide some com¬ 
ments on our numerical calculations. (*) Numerical cal¬ 
culations were performed for several values for the con¬ 
stituent quark mass M, from 350 MeV to 1000 MeV. 
The proper-time cut-off parameter A was not a free pa¬ 
rameter, but was determined so as to reproduce the pion 
decay constant /,r = 93 MeV[ p^. (ii) W e chose the initial 
functions Fo{p,z) = —tt -|- n^/f^ + ^/R, with R — 1.0, 
and 00 (p, z) = tan“^(p/z). We tried several forms of the 
initial functions and confirmed that the final result was 
independent of these choices, (m) Diagonalization of the 
Dirac hamiltonian was done following the method of Ka- 
hana and Ripka[p^. They used a discretized plane wave 
basis in a spherical box with radius D, which were defined 
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FIG. 2: Contour plot of the self-consistent profile function FIG. 3: Contour plot of the self-consistent profile function 
F{p,z), with M = 400 MeV. 0(p,z), with M — 400 MeV. 


by their grand spin and the parity. As previously stated, 
since our Dirac hamiltonian had axially symmetric prop¬ 
erty, the grand spin K was no longer good quantum num¬ 
ber of the eigenstates. In our case, since the third com¬ 
ponent of the grand spin and the parity tt = ± were 
only the good quantum number of the states, then we 
modified the Kahana-Ripka basis into those which were 
defined by the K 3 and the parity tt = ±. The new basis 
are the eigenstates of the free hamiltonian in a cylindri¬ 
cal box with height 2 x and radius Dp. These basis 
enable us to diagonalize our Dirac hamiltonian. 

In Figs. 2-3, we present results for the profile func¬ 
tions F{p,z), 0(p,z) with M = 400 MeV. In Figs. 4-5 
we display the baryon number density. Fig. Hshows the 
contribution from the valence quark and FigTB from the 
sea quark. From Figs. 4-5 it is found that the baryon 
number density has a toroidal shape. This result is con¬ 
sistent with other chiral invariant models using axially 
symmetric meson fields, such as the Skyrme model flltl 
and a naive quark meson model which involved six va¬ 
lence quarks and a pion cloud Q. The classical soliton 
energies corresponding to various values of M are given 
in Table ||. As M increases, the valence quark contribu¬ 
tion rapidly decreases, while that of the sea quark grows 
rapidly. Around M ^ 650 MeV, the valence level crosses 
zero energy and dives into the negative-energy region. 
For M > 650 MeV, the systems are dominated by the 


sea quark contribution. As for the total energy, which is 
sum of the valence quark and the sea quark contribution, 
there is essentially no noticeable change as M increases. 
This is a characteristic behaviour of our solution. The 
increase in the total energy value for large M is perhaps 
due to the omission of the higher wave numbers from our 
basis. 

Here, one could consider the classical “binding energy” 
, which is given by 


Dbound — Dstatic{U] ‘^^B — l.hedgehog- 


The classical nucleon energy with the hedgehog ansatz 
Mb=i hedgehog has been calculated by various authors 
p^ , pi] , [T^ , in which the values are chosen around 1200 
MeV. Thus, the classical binding energyEtound is about 
70 MeV at M = 400MeV. This is close to the Skyrme 
model result! 19 , and superior to the quark meson model 
value l|2^ which is 219 MeV. 


The mean radius of the toroid for the quark distribu- 
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FIG. 4: Baryon number density from valence quark contribu- FIG. 5: Baryon number density from sea quark contribution, 

tion, with M = 400 MeV. with M = 400 MeV. 


tion is estimated by 


(p)„ = 2^0 J Pdpdzd(pp(j)l{p,(p,z)(l)o{p,(p,z) ,{26) 

ip)o = 


X 


J pdpdzdpp(j)lip, <p,z)cl)^{p,ip,z), 


(27) 


(p) = {p)v + {p)o- 


(28) 


These values are also given in Table || and show a rapid 
decrease with increasing M. This is reasonable, because 
here M is regarded as the coupling-constant between the 
quark and pion, so the larger M means a stronger quark- 
pion interaction. The stronger interaction may produce 
more compact solitons. At M = 400 MeV, which may be 
a suitable choice, we obtained {p) = 0.672 fm which is in 
qualitative agreement with the Skyrme model value 0.78 
fm@. 

In summary, we have obtained the axially symmetric 
B — 2 soliton solution of the SU(2) yQSM. The solution 
was obtained in a self-consistent manner. The results 
are in qualitative agreement with those from the Skyrme 
model and other quark meson models. This suggests that 
these features are independent of the particular choice of 
chirally invariant model. Individualities of each model 
will become clearer after thorough investigations for var¬ 


TABLE I: Glassical 
of toroid. 

mass spectrum (in 

MeV) and 

mean radius 

M 


E field 

Estatic 

(P) [fm] 

350 

1134 

1189 

2323 

0.705 

400 

925 

1397 

2322 

0.672 

450 

765 

1569 

2334 

0.629 

500 

558 

1760 

2318 

0.600 

600 

184 

2153 

2337 

0.549 

700 


2384 

2384 

0.508 

800 


2466 

2466 

0.482 

900 


2589 

2589 

0.462 

1000 


2763 

2763 

0.447 


ious physical observables . The most striking differ¬ 
ence between our yQSM and the Skyrme model is the 
existence of quark degrees of freedom. From consider¬ 
ation of the single quark energy level, we confirm that 
the minimum energy configuration of i? = 2 is axially 
symmetric, while in the Skyrme model it is a conjecture. 
Since the xQSM includes the valence and the sea quark 
degrees of freedom, we can give theoretical support for 
nuclear medium effects such as the EMC effect in deep 
inelastic scattering experiments. 

The solutions obtained here were classical ones which 
have no definite spin, isospin quantum numbers corre- 
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spending to physical particles. Therefore the solutions 
should be quantized by projecting onto good spin, 
isospin states in order to estimate the energies, the mean 
square radius, and other static properties of the physical 
B = 2 system. The quantization of our solution using 
the well known cranking procedure in SU(2) is now in 
progress. 

One of us (N.S.) is very grateful to Dr. S. Akiyama 
for many valuable discussions and comments from the 
beginning of this work. 
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